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Abstract 


Optical  depth  measurements  of  transmissive  cirrus  clouds  were  made  using  coincident  lidar 
and  satellite  data  to  improve  our  interpretation  of  recent  satellite  cloud  climatologies.  These 
climatologies  differ  in  the  way  they  detect  transmissive  clouds  because  some  use  solar  reflectance 
data  (ISCCP)  while  others  use  multi-spectral  infrared  data  (C02  Slicing).  To  relate  these 
climatologies  and  estimate  the  impact  of  transmissive  clouds  on  the  earth’s  heat  budget,  a 
relationship  between  visible  and  infrared  radiation  properties  has  to  be  used.  We  examined  the 
popular  assumption  that  the  ratio  of  the  visible  to  infrared  optical  depths  should  be  2.0  because  the 
visible  extinction  cross  section  is  twice  the  infrared  absorption  cross  section  when  cloud  particles 
are  large  compared  to  the  wavelength. 

Our  lidar-satellite  data  combination  confirmed  this  relationship  for  optically  thin  cirrus 
(Tvis<  1)  that  were  uniform  both  vertically  and  horizontally.  Deviations  from  the  2.0  ratio  appeared 
to  come  from  differences  in  sampling  horizontally  variant  clouds.  The  lidar  observed  a  small  area 

A 

.  (~i  m  diameter)  with  the  only  spatial  sampling  being  from  wind  advection  of  the  cloud  over  the 
sensor  while  the  satellite  sensors  sampled  larger  areas  1-20  km  in  diameter  in  both  the  along  wind 
and  cross  wind  directions. 

Most  cirrus  were  physically  thick  (  >  1km)  and  often  composed  of  multiple  layers.  Super 
cooled  liquid  water  layers  embedded  in  ice  clouds  were  found  in  32%  of  the  cirrus  examined. 

For  thick  cirrus  and  multi-layered  cirrus,  a  direct  optical  depth  comparison  was  not 
possible  because  of  the  single  layer  model  used  in  the  satellite  analysis.  However,  we  simulated 
what  the  satellite  infrared  analysis  (C02  Slicing)  would  produce  in  complicated  cirrus  structures  of 
multiple  levels  including  a  super  cooled  liquid  water  layer.  The  simulation  used  the  detailed 
measurements  of  vertical  extinction  from  the  lidar  for  deriving  the  radiances  that  the  satellite 
would  receive  on  each  of  its  infrared  channels.  A  single  cloud  height  solution  was  calculated 
following  the  C02  Slicing  method.  An  infrared  optical  depth  was  also  calculated  for  this  single 
cloud  level  using  the  satellite  measured  radiances  in  the  10.8  /im  window  channel.  The  calculated 
cloud  level  was  close  to  the  altitude  of  the  densest  cloud  layer  found  by  the  lidar.  The  infrared 
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optical  depth  (10.8pm)  was  =0.5  of  the  visible  optical  depth  (0.53  pm)  integrated  over  all  the  cloud 
layers  which  agreed  with  theory.  This  indicates  that  the  visible  and  infrared  satellite  cloud 
climatologies  derived  from  can  be  related  for  long  term  studies. 
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1.  INTRODUCTION 


Predictions  of  global  warming  from  the  increase  of  C02  in  the  atmosphere  have  lacked 
credibility  because  of  the  need  to  assess  the  contribution  of  clouds  to  the  radiation  budget 
(Paltridge,  1980).  Cirrus  clouds,  in  particular,  have  the  ability  to  either  offset  the  C02  warming  by 
increased  solar  reflection  or  increase  the  warming  by  the  capture  of  terrestrial  radiation.  Thin 
cirrus  in  the  upper  troposphere  will  add  to  tropospheric  warming  because  it  allows  solar  insolation 
while  capturing  par,  of  the  terrestrial  radiation.  Thick  cirrus  and  multiple  layer  cloud  conditions 
where  water  clouds  are  present  under  the  cirrus  produce  highly  reflective  cloud  combinations  tha, 
can  lower  albedo  and  solar  insolation.  Predictions  of  climate  changes  ultimately  must  include  the 
radiative  effects  of  these  clouds.  However,  Stephens  et  al.  (1990)  point  ou,  that  the  effect  of 
clouds  on  the  earth’s  radiation  budget  is  not  understood  because  of  the  lack  of  data. 

G16ba!  cloud  and  radiation  data  sets  are  presently  being  collected  from  satellites.  The 
International  Satellite  Cloud  Climatology  Project  (ISCCP,  Schiffer  and  Rossow,  1983,  and  Rossow 
and  Schiffer,  1991)  has  measured  cloud  cover  and  radiation  for  1 1  years  using  the  operational 
weather  satellites.  Wylie  e.  al.  (1994)  have  studied  citrus  clouds  from  four  years  of  polar  orbiting 

NOAA  satellite  data,  also  Woodbury  and  McCormick  (1983)  have  made  a  similar  analysis  of 
global  cirrus  from  the  SAGE  H  satellite  data. 

Inter-annual  changes  in  cloud  cover  have  been  found.  The  Wylie  et  al.  (1994)  study 
reported  an  increase  in  cirrus  with  the  recent  El  Nino-Southem  Oscillation  (ENSO)  tha,  began  in 
1991.  A  similar  change  in  cloud  cover  in  the  1982-83  ENSO  was  also  found  by  Weare  (1992)  using 
the  NIMBUS  7  data  set.  TTtese  studies  indicate  tha,  doud  cover  is  linked  to  climatic  fluctuations. 

Estimates  have  been  made  of  the  effect  of  lower  tropospheric  clouds  on  the  earth's  hea, 
budge,  by  Tselioudis  e,  al.  (1993)  using  ISCCP  data,  bu,  similar  calculations  for  upper 
tropospheric  clouds  have  not  been  attempted  because  of  the  complications  of  estimating  the 
effects  of  both  the  visible  and  the  infrared  radiative  processes  together.  Boundary  layer  and  lower 
tropospheric  douds  affect  solar  insolation  bu,  have  little  effect  on  infrared  cooling,  while  cirrus 


clouds  affect  both  radiative  processes.  Satellite  observations  have  usually  been  restricted  to  either 
the  visible  or  the  infrared.  The  satellite  data  used  by  Wylie  et  al.  (1994)  and  Weare  (1992)  were 
from  longwave  infrared  sensors,  therefore  the  effect  on  solar  radiation  will  have  to  be  estimated 
smce  it  was  not  measured.  The  ISCCP  differs  from  these  two  studies  in  that  it  employs  visible 
sensors  for  detecting  transmissive  cirrus  clouds.  The  IR  contribution  must  be  estimated  from  the 
visible  reflectance  measurements.  The  11  year  record  of  the  ISCCP  (Schiffer  and  Rossow,  1983) 
will  allow  studies  of  fluctuations  of  cloud  cover  with  many  short  term  climate  changes  including 
three  ENSO’s.  The  SAGE  H  data  set  can  also  be  used  for  deriving  visible  optical  depths  (Law,  et 
al.  1994a  and  b).  These  satellite  data  sets  are  complimentary  where  they  overlap,  but  the  overlap 
is  small,  about  one  year  for  each  intersection.  To  study  interannual  changes,  trends,  and  longer 
term  fluctuations,  each  of  these  data  sets  will  have  to  be  used  separately.  Visible  to  infrared 

radiance  relationships  will  have  to  be  used  to  fully  estimate  the  effect  of  cirrus  changes  on  heat 
'*  budget  parameters. 

Our  knowledge  about  the  relationship  of  the  solar  reflection  to  the  infrared  attenuation  in 
cirrus  clouds  is  limited.  For  particles  which  are  large  compared  to  the  wavelength,  we  expect  the 
visible  scattering  cross  section  to  be  twice  the  infrared  absorption  cross  section  when  there  is  little 
visible  absorption  and  strong  infrared  absorption  by  ice  and  liquid  water.  This  implies  that  the 
ratio  of  the  visible  to  the  infrared  optical  depth  should  be  about  2.0.  Modelling  of  the  radiative 
properties  of  ice  crystals  by  Minnis  et  al.  (1993)  and  Mitchell  and  Amott  (1994)  support  this 

assumption.  However,  variances  in  cloud  particle  size  and  shape  could  vary  this  ratio  from  1.8  to 
4.0. 

In  situ  measurements  of  cirrus  microphysical  properties  have  found  extremely  large 

variances  in  number  density,  ice-water  content  and  particle  size  (Dowling  and  Radke,  1990).  The 

number  density  of  cloud  particles  varies  by  8  orders  of  magnitude,  from  10*4  L*1  to  10 +4  L*1.  The 

ice-water  content  varies  over  10  g  nr3  and  the  particle  size  can  vary  from  1  to  8000  pm.  With  this 

wide  range  of  possible  microphysical  variances,  large  variances  in  the  radiative  properties  are 
expected. 


Few  measurements  of  cirrus  have  been  made  when  data  from  both  the  visible  and  IR  were 
collected  simultaneously.  Plat.  (.979)  used  a  vertically  pointing  ,R  radiometer  with  a  visible  Udar 
to  study  a  limited  number  of  thin  cirrus  cases.  Minnis  e«  ai.  (1990)  studied  one  case  from  FIRE-f 
(the  First  ISCCP  Regional  Experiment  in  October  1986)  using  infrared  data  from  a  satellite  and 
visible  data  from  a  lidar.  Plat,  (1979)  found  a  visible  to  IR  optical  depth  ratio  of  2.0  while  Mimes’ 
case  was  slightly  higher  at  2.1.  These  cases  were  mostly  thin  single  layer  cirrus. 

To  use  the  global  cloud  climatologies,  we  need  to  understand  what  has  been  measured  and 
how  their  representation  of  clouds  depicts  the  radiative  transfer  in  both  the  visible  and  the 
infrared.  The  retrieval  of  cloud  parameters  in  the  satellite  climatologies  has  been  based  on  the 
assumption  that  the  clouds  are  a  single  thin  layer  and  that  they  cover  the  Field  of  View  (FOV)  of 
the  satellite  sensor  for  each  pixel  in  the  image.  ISCCP  and  Wylie  et  al.  (1994)  only  solved  for  one 
cloud  level  in  each  pixel  (FOV).  A  comparison  of  satellite  cloud  retrievals  to  ground  based 
Rational  Weather  Semce  obsetvations  over  the  continental  United  States  found  that  low  and 
middle  tropospheric  clouds  were  present  in  52%  of  the  cases  where  the  satellite  cloud  detection 
system  found  an  upper  tropospheric  cloud  (Menzel  et  al,  1992).  Similar  statistics  on  the  frequency 
of  multiple  cloud  layers  over  oceans  are  shown  in  Hahn  et  al.  (1982).  Lidar  data  collected  at  the 
University  of  Wisconsin-Madison  (Figure  1)  found  that  cirrus  are  often  geometrically  thick 
reaching  depths  up  to  5  km.  Thus,  we  expect  that  cirrus  are  more  complicated  than  the  assumed 
thrn  single  layer  model  used  in  the  cloud  retrieval  algorithms. 

Because  of  the  uncertainty  of  wha,  satellite  cioud  retrievals  represent  and  the  need  to 
relate  the  different  climatologies  so  .ha,  longer  cloud  records  can  be  compiled,  we  collected  a  data 
se,  of  high  quality  lidar  observations  with  coincident  satellite  data.  The  purpose  of  this  study  is 
twofold:  1)  to  asses  the  applicability  of  the  ’sin^e  thin  layer-  assumption  used  in  satellite  cloud 
retrievals,  and  2)  to  asses  the  validity  of  the  2.0  visible  to  htfrared  optical  depth  ratio  so  the 
satellite  cloud  climatologies  can  be  combined  for  a  more  complete  assessment  of  radiative  transfer 


processes. 


2.  MEASUREMENT  TECHNIQUE 

™*  S,Udy  WaS  “  1989  usi"8  the  unique  measure*  systems  a.  ,h.  University  of 

Wisconsin-Madison  (Ackerman,  e.  al„  1993).  Visible  optical  depths  of  cirrus  were  obtained  from 
the  High  Spectral  Resolution  Lidar  (HSRL)  while  mfrared  optical  depths  were  obtained  from 
satellite  data  so  that  the  information  collected  could  be  related  i 
et  al.  (1994)  and  ISCCP. 

2jh.  The  HSRL 


to  the  cirrus  climatologies  of  Wylie 


In  order  to  make  reliable  measurements  of  the  extinction  profile,  the  HSRL  measures  two 
signals  from  the  0.53  am  transmitted  pulse  which  can  be  processed  to  yield  separate  lidar  returns 
for  aerosol  and  molecular  scattering.  This  separation  is  possible  because  the  spectrum  of  the 
molecular  lidar  return  is  Doppler  broadened  by  the  thermal  motion  of  the  molecules  while  the 
slow  moving  aerosol  panicles  generate  neghgible  spectral  broadening.  The  molecular  scattering 
JOSS  section  is  a  function  of  molecular  density  and  can  be  calculated  from  Rayleigh  scattering 
theory  and  an  independently  measured  temperature  profile.  The  HSRL  then  uses  molecular 
scattering  as  a  calibration  target  which  is  available  at  each  point  in  the  lidar  return. 

Two  lidar  returns  are  derived;  these  are  given  by  the  following  equations  for  the  signal 
power  received  from  molecular  scattering,  Pm(r),  and  aerosol  scattering,  Pa(r): 


P«(r)  ~  Eo  2 ?!  M r)  ~g^~exP(m2r(r)) 


P»  -  Eo  Pa( 0  exp(-2r(r)) 

Eo  =  Laser  pulse  energy,  J. 
c  =  Speed  of  light,  ms-1. 

Ar  =  Collecting  area  of  the  receiver,  m2. 
r  -  Range  to  the  scattering  volume,  m. 
ia  2  ^ert?S0.1  scattering  cross  section  per  unit  volume,  m-\ 

P(ff>r )  ttermg  Cross  sectlon  Per  unit  volume  for  air  molecules,  m-l. 

=  Backscatter  phase  function,  sr1. 


(1) 

(2) 


-1 


4?r 


-  Optical  depth. 


7 


Equation  ,  contains  two  unknowns:  the  optical  depth  between  the  lidar  and  the  scatter 

volume,  r  (r),  and  the  molecular  scattering  cross  section  a  M  Mol  l 
bv  sectl°n./«r).  Molecular  scattering  is  described 

ay  etgh  scattermg  equation  and  is  directly  proportional  ,o  atmospheric  density,  ,(r). 

^m(r)  =  Cp(r) 

Thus  if  we  define  a  range-squared,  energy  corrected  molecular  lidar  return: 

Sm(r)  =  Pm(r)r2/E0  ’  ' 

**Pm(r)]  =  constant  +  ln[p( r)]  .  2t(t) 

n.e  optical  depth  between  any  two  ranges,  r  and  r„,  is  given  by: 

r(r)',(r»)  =  a5feWr)/,M  -  0.51n[Sm(r)/Sm(ro)] 

The  average  extinction  cross  section  in  a  layer  between  r  amt  r  k 

J >  (r(r)  -  r(r0)]/[r-rj  °  C°mpU,ed  from  3: 

W 


(1994). 


More  details  can  be  found  in  Grund  and  Horan, a  (1991)  and  Piironen  and  Horan, a 


2.b*  Satellite  Infrared 


The  satellite  cloud  analysis  of  Wylie  et  al.  (1994)  assumes  that  the  radiative  effects  of  the 

clouds  occur  at  only  one  level  in  the  troposphere  This  all  •  , 

P  P  This  allows  a  simple  calculation  of  cloud 
eminence  a,  10.8  pm,  the  longwave  window  channel,  from  the  radiances 
in  the  cloudy  pixel  (R^),  a  clear  pixe|  ( ^  an(J  ^ 

present  at  the  pixel  (Rc]d). 


measured  by  the  satellite 
radiance  expected  if  an  opaque  cloud  were 


'a  -  tRd.-RaJ/tRa.-R.J 
Rut  was  measured  from  the  pixel  that  covered  the  HSRL  while  is  an  estimate  of  the  radiance 

derivedf  J ^ ^ meaS“ed “ ‘he d°Ud "°‘ PreSe"‘'  ^ radiance «“■»«« was 
derived  from  the  nearest  clear  pixel  by  assuming  that  the  background  under  the  cloud  was  the 

same  as  at  the  location  where  Rse  was  taken.  Rua  was  calculated  from  the  temperature  of  the  mid 
cloud  level  using  rawinsonde  data  from  Green  Bay,  WI  which 


(5) 


was  the  nearest  sounding  to  the 
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HSRL.  To  simplify  the  calculation,  the  black  body  temperature  to  the  fourth  power  (T4)  was  used 
in  place  of  R.  The  satellite  measurements  are  calibrated  to  both  radiance  (R)  and  blackbody 
temperature  (T).  Equation  5  was  expressed  as: 

fir  =  [T4clr-,T»sat]/[nfr-T4cld]  (6, 

where: 

Tcir  =  blackbody  brightness  temperature  for  a  clear  pixel 
Tsat  =  blackbody  brightness  temperature  for  a  cloudy  pixel 

Tdd  =  the  air  temperature  at  the  mid-cloud  height. 

The  mid-cloud  altitude  was  determined  from  the  appearance  of  the  cloud  on  the  HSRL 
backscatter  vertical  profile.  This  approximated  the  cloud  level  that  the  HIRS  analysis  reports  for 
transmissive  clouds  (Wylie  and  Menzel,  1989). 

IR  optical  depths  were  calculated  from  the  emittances  assuming  no  scattering  at  the 
window  channel  IR  wavelength  (10.8  /im)  and  that  transmission  was  the  compliment  of  emission: 

;  *  ‘  =  ’Mi-fir) 

The  infrared  optical  depths  were  also  corrected  for  the  angle  of  the  satellite  scan  through  the 
atmosphere  using  the  cosine  of  the  zenith  angle. 


3.  DATA 

HSRL  data  were  examined  for  21  separate  days  covering  all  four  seasons  (see  Table  1). 
NOAA  12  overpasses  of  Madison,  WI  were  available  when  most  of  the  HSRL  data  were  taken. 
However,  on  six  days,  GOES  7  data  were  compared  to  the  HSRL  because  of  the  absence  of 
NOAA  12  data.  Comparisons  were  restricted  to  cases  with  visible  optical  depths  <  3  when  the 
HSRL  probed  through  the  entire  cloud.  Total  cloud  penetration  was  determined  by  the  presence 
of  molecular  backscatter  signal  beyond  the  cloud  top  (Equation  1). 


4.  RESULTS 

4a  Visible  and  IR  Optical  Depths 
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To  compare  visible  and  infrared  radiative  properties,  satellite  infrared  optical  depths  from 
Equation  7  were  plotted  with  the  HSRL  visible  optical  depths  from  Equation  3  (see  Figure  2). 
The  21  cases  covered  a  variety  of  transmissive  cirrus  situations.  The  clouds  were  separated  into 
three  classes:  1)  geometrically  thin,  covering  <  2.5  km,  2)  geometrically  thick,  covering  >  2.5  km, 
and  3)  multi-layered  clouds.  Many  of  the  multiple  level  cloud  cases  included  a  thin  layer  of  super 
cooled  liquid  water  (distinguished  by  low  values  of  depolarization  measured  by  the  HSRL)  and 
often  included  ice  crystals  precipitating  from  the  layer.  Liquid  water  layers  embedded  in  cirrus 
were  found  on  32%  of  the  days  that  the  HSRL  collected  cirrus  data. 

Error  estimates  for  both  the  HSRL  and  satellite  optical  depth  measurements  are  also 
included  in  Figure  2.  For  the  visible  optical  depth  (horizontal  bars)  the  error  was  computed  from 
an  error  analysis  including  photon  counting  statistics,  error  in  the  observed  density  profile,  and 
error  in  system  calibration.  For  the  satellite  infrared  optical  depths  (vertical  bars),  the  error  was 
^sfimated  by  assuming  a  2  K  error  in  each  variable  in  Equation  6.  This  is  our  estimate  of  the 
maximum  uncertainty  in  determining  the  clear  radiance  (T»clr)  and  the  cloud  level  radiance  (T>dd). 
The  largest  source  of  clear  radiance  error  was  spatial  variances  or  gradients  in  the  land  surface  or 
low  level  aerosol  density  between  the  locations  of  the  HSRL  and  where  Tclr  was  obtained.  The 
largest  source  of  Tcld  error  was  the  air  temperature  gradient  between  the  location  of  the  HSRL 
and  the  rawinsonding.  The  satellite  radiometric  measurements  are  calibrated  to  0.8  K  which  is  a 
far  better  accuracy  than  the  other  error  sources.  The  location  of  the  HSRL  in  the  satellite  images 

was  also  determined  to  an  accuracy  of  <  1  pixel  using  navigation  adjustments  from  rivers  and  lake 
shorelines. 

4.b.  Single  Layer  Cirrus 

An  example  of  a  single  layer  cirrus  cloud  is  shown  in  Figure  3  from  26  October  1993.  The 
cirrus  cloud  was  found  between  9  and  11  km  and  was  fairly  uniform  in  the  vertical  dimension  and 
in  time  as  seen  in  the  lidar  time  section.  The  NOAA  satellite  image  (Figure  4)  also  showed  very 
uniform  conditions  around  the  HSRL  (black  dot  on  the  image).  Brightness  variations  around  the 
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HSRL  Were  hiSHy  "  *“  a  °f°"*  > «  »  blackbody  temperature.  ^  cim$ 

isible  optical  depth  was  0.5  at  the  time  of  the  NOAA  ovetpass  (0:30  UT).  The  satellite  infrared 

opnca,  depth  was  found  to  be  0.28,  yieldmg  a  Wrt  ratio  of close  to  the  expected  2.0 

Low  level  aerosols  were  also  present  from  0-3.2  km  as  indicated  by  the  HSRL  (Fig.  3) 

These  aerosols  were  probably  not  included  in  the  satellite  optical  depth  measurement.  Tc(r  in 

Canon  5  was  taken  from  an  area  ,50  km  to  the  southeast  in  the  sateUite  image  where  similar 
low  level  aerosols  should  have  been  present. 


4.c.  The  Effects  of  Multiple  Cloud  Layers 

Multiple  layers  of  clouds  have  been  vety  common  during  this  study.  Five  cases  of  multiple 
layers  were  examined  in  this  study.  One  example  is  shown  in  Figure  5  for  11  November  ,993  when 

VT“S  ”S  dUrin8  ‘he  HSRL  °bSerVa'“n-  one  layer  of  cirrus  from  7  to 

PreSent* tben  at  3  ,econd  layer  appeared.  Depolarization  measurements  by 

=  mdicated  that  the  mtense  backscatter  levels  (solid  white)  were  liquid  water  particles 

"  6  “  °f  ‘he  Cl°Ud  fce-  Visibk  <**•«  depths  increased  from  0.4-0.6  to  1.0  a,  the 

NOAA  ovetpass  at  1:30  UT.  After  1:40  UT  the  visible  optical  depth  exceeded  3.0,  the  maximum 

measurable  by  the  HSRL.  The  NOAA  12 IR  image  in  Figure  6  also  shows  this  change  as  dense 
C  OU  s  moved  in  from  the  west.  Notice  that  the  light  cirrus  to  the  east  (down  wind)  has  faint 
northwest- southeast  stria, ions  while  the  denser  area  to  the  west  (up  wind)  has  onhagonal 
sou, hwes, -northeast  stria, ions.  Both  patterns  are  skewed  to  the  westerly  wind. 

"Hie  water  layer  dominated  the  optica,  depth  of  the  combined  cirrus  layers.  Theupperice 
cloud  maintained  a  nearly  constant  optical  depth  of  0.4-0,  until  01:40  UT  when  the  HSRL  could 
no  longer  penetrate  through  the  water  layer. 

The  satellite  measured  blackbody  radiance  decreased  from  261 K  before  the  lower  layer 

n  ^to"K“‘"HsRLk^ 

roppe  to  243  K  at  30  km  west  of  the  HSRL  in  the  image. 
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Simulations  of  the  HIRS  radiance  and  cloud  height  measurements  were  made  for  three 
situations:  1)  the  single  layer  case,  2)  the  two  layers  at  the  NOAA  12  over  pass  at  1:30  UT,  and  3) 
the  point  when  the  lower  water  layer  first  obscured  the  HSRL  at  1:40  UT.  The  second  and  third 
situations  represent  problems  to  the  basic  assumption  used  in  satellite  cloud  retrievals  -  the  single 
level  description  of  the  cloud.  The  expected  response  of  the  multi-spectral  HIRS  cloud  retrieval 
algorithm  used  by  Wylie  et  al.  (1994)  to  multi-layered  clouds  is  that  it  would  place  the  single  cloud 
height  between  the  two  layers  (Wielicki  and  Coakley,  1981  and  Menzel  et  al.,  1992).  The  single 
level  cloud  height  is  calculated  from  the  mean  altitude  of  infrared  radiative  divergence,  similar  to  a 
center  of  mass  concept  in  kinematics.  However,  when  one  cloud  layer  has  more  optical  depth  than 
the  other,  the  multi-spectral  HIRS  cloud  height  solution  is  expected  to  move  toward  the  denser 
layer.  An  exception  to  this  trend  occurs  when  the  top  layer  is  nearly  opaque  and  obscures  the  view 
of  the  lower  layer,  forcing  the  HIRS  solution  into  the  top  layer. 

1  '  '  AttemPts  have  been  made  to  extract  more  than  one  cloud  level  from  the  multi-spectral 
HIRS  data.  However,  the  solution  of  the  radiative  transfer  model  for  the  HIRS  sensors  requires 
the  integral  of  the  radiation  over  a  broad  portion  of  the  troposphere.  The  cloud  height  solution 
has  been  mis-labelled  the  C02  Slicing  Method  from  the  desire  that  each  channel  would  see  a 
different  "slice"  of  the  atmosphere;  however,  altitude  sensitivities  of  the  channels  are  broad  and 
overlap.  Each  HIRS  channel  receives  radiation  from  3-5  km  of  the  troposphere.  To  solve  for 
more  than  one  cloud  level,  the  vertical  profile  of  the  clouds  must  be  determined.  Since  each  HIRS 
sensor  sees  the  integral  of  the  radiation  over  a  large  depth  of  the  troposphere,  a  wide  variety  of 

cloud  structures  could  be  present  from  which  no  unique  solution  can  be  made.  The  only  tractable 
solution  for  the  HIRS  data  is  one  cloud  level. 

The  response  of  the  multi-spectral  HIRS  algorithm  to  the  cloud  structures  observed  by  the 
HSRL  was  simulated  by  calculating  what  the  HIRS  channels  should  have  seen  given  the 
temperature  structure  and  the  infrared  emittance  and  transmittance  of  each  cloud  layer.  HIRS 
data  were  available  for  all  cases  where  the  NOAA  satellites  were  used.  The  HIRS  data  were  not 
used  because  the  HIRS  data  were  not  directly  over  the  HSRL.  The  HIRS  pixels  (FOV)  do  not 
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overlap  and  contain  gaps  between  pbtels  tha,  are  about  the  same  dimensions  as  the  FOVs.  For 

urate  comparisons  to  the  HSRL,  radiative  data  must  come  from  the  location  of  the  HSRL 

patral  vanance  of  the  clouds  can  cause  large  errors  These  errors  are  discussed  in  derail  in 
section  4d. 

Tbe  radiance  received  by  each  HIRS  sensor  was  esrimared  using  .  simple  mode,  of  rhree 
cloud  layers. 

**0 


-  JpiB[1.T(p)]{dSi(p)/dp}dp  +  Sj^pi)  e  1B[i,T(pi)] 

+  (1*ei)[  Jp2B[i.T(p)]{dSi(p)/dp}dp  +  Si(p2)e2B(i,T(p2)]] 

+  (1  el)(i-e2)[  jp3B[1>T(p)]{dSi(p)/dp}dp  +  Si(p3)e3B(i,T(p3)]] 

+  (i  £i)(l-e2)(l-e3)[  i^sfcB[i,T(p)]{dSi(p)/dp}dp  +  S^sfcJBftTYsfc)]] 

Where:  J 

s$fL:  zdouezrr ^ 

€*  =  the  emittance  for  cloud  layer  x.  1  p  through  the  atmosphere, 

The  emittancc  (<)  for  ^  c)oud  )ayer  ws  approximated  (rom  ^  hsrl  meas^nK  ^ 

Vm  0P'iCai  deP'h  °f  ^ layer'  ^  «M.  was  assumed  ,o  be  one  ha,  of  rhe 

^  dCP,h  f°r  “*  d°Ud  ^  -  emirtance  W  was  derived  from  Equa,io„  7 

TTre  W~d  radiances  for  rhe  HIRS  channels  were  used  in  rhe  algorirhm  desorbed 
y  y  e  et  al.  (1994)  for  determining  a  single  cloud  level.  This  cloud  level  was  then  used  for 
etermining  TcJd  from  the  rawinsonde  data.  With  a  new  T  anH 

wun  a  new  ic]d  and  the  measured  T  andT 

Equations  6  and  7  were  re-applied  to  give  an  estimate  of  wha,  the  HIRS  algorithm  would  L 

enved  for  the  cloud  altrtude  and  infrared  optical  depth  (rt).  Note  that  this  calculation  used  a 

com  mat, on  of  satellite  measured  radiances,  Tsat  and  Tclr  and  a  new  radiance  for  the  approximate 
single  cloud  level,  Tcld. 


(8) 
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The  first  situation  on  11  November,  was  a  single  cirrus  layer  with  a  mean  height  of  8.5  km. 
The  infrared  optical  depth  calculated  from  Equations  5  and  6  was  0.25.  With  the  HSRL  visible 
optical  depth  measurement  of  0.5,  the  T^/ra  ratio  of  2.0  agreed  with  theory. 

The  HIRS  radiances  for  the  second  situation  at  1:30  UTwere  calculated  from  Equation  8 
using  two  cloud  levels  at  5  and  8.5  km.  The  visible  optical  depths  of  the  two  layers  measured  by  the 
HSRL  were  0.9  (lower)  and  0.6  (upper).  Emittances  of  0.4  for  the  lower  layer  and  0.3  for  the 
upper  layer  were  calculated  from  Equation  7  by  assuming  that  the  infrared  optical  depths  were  one 
half  of  the  visible  HSRL  measurements.  The  single  cloud  layer  calculated  from  the  Wylie  et  al. 
(1994)  algorithm  was  6.6  km  which  was  between  the  two  cloud  layers.  The  infrared  optical  depth 
from  this  level  and  the  AVHRR  T^  and  Tclr  measurements  was  0.6.  Using  the  HSRL  measured 
visible  optical  depth  of  1.0,  Tvv/Ta  was  1.7,  slightly  lower  than  theory. 

In  the  third  situation  at  1:40  UT  the  total  optical  depth  exceeded  the  HSRL  measurement 
'1  capability.  However,  Equation  8  was  employed  to  estimate  HIRS  radiances  by  assuming  that  the 
upper  cloud  layer  remained  at  the  same  optical  depth  as  measured  at  1:30  UT  and  the  lower  cloud 
layer  was  at  the  level  of  HSRL  saturation,  3.0  visible  or  1.5  infrared.  The  single  cloud  level 
solution  was  5.9  km,  0.7  km  closer  to  the  liquid  water  layer  than  the  calculation  from  the  second 
situation. 

This  is  an  example  of  how  the  HIRS  single  cloud  level  solution  used  by  Wylie  et  al.  (1994) 
responds  to  multiple  level  clouds.  The  technique  is  constrained  to  reporting  only  one  cloud  layer. 
For  one  layer  of  transmissive  cirrus,  it  places  the  cloud  level  near  the  center  of  the  layer  as 
predicted  by  Wielicki  and  Coakley  (1981).  Similar  results  for  a  previous  comparison  with  lidar 
data  were  reported  in  Wylie  and  Menzel  (1989).  For  two  layers,  the  HIRS  solution  was  between 
the  layers,  moving  toward  the  denser  lower  layer  as  it  became  the  dominant  layer. 

We  acknowledge  that  the  estimation  of  HIRS  infrared  radiances  and  the  calculation  of 
infrared  optical  depths  for  the  two  layered  cloud  appears  to  be  a  circumlocution  because  HSRL 
visible  optical  depth  measurements  were  used.  However,  the  infrared  optical  depth  calculation 
(Equations  6  and  7)  involves  more  than  the  cloud  height  and  the  radiance  (T4^  from  this  cloud 
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height  level  using  .he  temperature  sounding.  The  other  terms  in  Equation  6,  Tsal  and  Tclp  came 
from  satellite  measurements  independently  of  the  HSRL  measurements.  Also.  Equation  8  uses 
radiances  from  the  sounding  (the  integrals)  as  well  as  the  eminences  of  the  cloud  layers  .ha,  came 
from  the  HSRL  optical  depths.  To  test  the  sensitivity  of  Equation  8  to  the  method  of  transfotmmg 
HSRL  visible  optical  depth  measurements  to  infrared  eminences,  the  calculations  were  repeated 
using  Wm  ratios  of  1.0  and  4.0  on  the  HSRL  visible  optical  depth  measurements.  These  two 
tests  differed  from  the  original  relationship  by  a  factor  of  2.  These  tests  changed  the  single  cloud 
level  solution  by  ±0.6  km  and  the  derived  infrared  optical  depth  by  ±0.1.  This  is  smaller  than  the 
changes  caused  by  the  increase  in  density  of  the  lower  cloud  layer. 


4.d.  The  Effects  of  Vertical  and  Horizontal  Variability 

Most  of  die  deviations  from  the  2.0  visible  to  IR  optical  depth  ratio  in  Figure  2  came  from 
cirrus  clouds  with  high  horizontal  variations.  An  example  of  such  a  cloud  observed  with  the  HSRL 
•on  2  September,  93  is  shown  in  Figure  7.  The  NOAA  12  satellite  image  from  the  01:36  UT 
overpass  (Figure  8)  shows  a  small  and  long  cirrus  cloud  which  is  partially  over  the  HSRL  (dot  in 
•he  image).  This  cloud  moved  east  on  a  75»  bearing.  Only  part  of  it  crossed  the  HSRL.  Larger 
clouds  up  wind  crossed  the  HSRL  a,  01:55  UT.  This  is  an  example  of  how  spatial  variability 
affected  the  lidar-satellite  comparison.  The  visible  optical  depth  a,  the  time  of  the  satellite 
overpass  was  0.2  while  the  infrared  optical  depth  a,  the  presumed  HSRL  location  was  0.05  yielding 
a  Wm  ratio  of  4.0,  larger  than  predicted  from  theory.  The  satellite  radiance  was  taken  in  a 
gradient  region  where  any  error  in  location  could  raise  or  lower  its  value.  Also,  for  horizontally 
smaU  and  narrow  clouds,  the  larger  size  of  the  satellite  radiometer’s  FOV  may  have  averaged  more 
cloud  area  with  less  cloud  optical  depth  than  that  viewed  by  the  HSRL. 

Most  of  the  low  optical  depth  cases  shown  in  Figure  2  came  from  GOES  satellite  data 
which  has  a  larger  FOV  than  the  AVHRR  (10.8  km  diameter  GOES  vs.  the  1  km  AVHRR).  The 
GOES  tended  to  find  higher  optical  depths  than  predicted  by  the  lidar,  the  opposite  of  the  2 
September  case  previously  discussed.  Substantial  horizontal  variances  were  present  in  all  of  the 
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cases  with  visible  optical  depth  <  0.2.  We  suspect  that  the  satellite  pixels  used  to  determine  clear 
radiances  (T**)  may  have  had  cloud  contamination  and  thus  may  have  been  too  cold  which 
would  cause  higher  infrared  optical  depths  to  be  calculated. 

Another  example  of  a  variable  cloud  situation  was  found  on  15  April,  94,  see  Figures  9  and 
10.  ITte  HSRL  image  of  backscatter  cross  section  (Figure  9)  showed  large  variations  in  the  optical 
depth  with  time  while  the  NOAA  12  AVHRR  image  (Figure  10)  showed  cirrus  clouds  with  a 
variety  of  shapes.  An  elongated  cloud  partially  covered  the  HSRL  at  the  time  of  the  ....no. 
overpass  (00:49  UT)  and  moved  toward  the  northeast  on  a  80°  bearing  There  were  at  least  two 
layers  of  cirrus  as  indicated  in  Figure  9.  The  lower  layer  was  more  dense  but  highly  variable  in 
nature  and  contained  some  super  cooled  water.  The  upper  layer  was  very  thin  and  composed 
entirely  of  ice.  He  visible  optical  depth  was  0.75  while  the  infrared  was  1.07.  This  combination 
produced  ata/t,  optical  depth  ratio  of  0.7,  far  less  than  the  theoretical  2.0. 

?  ■  ■  T0lla*,iCaliy  eS'ima,e  ,he  OTOrs  °fFOV  ske  differences  and  possible  miss-registration 
of  the  satellite  image,  data  from  the  Volume  Imaging  Lidar  (VIL)  were  examined.  This  is  a  1064 

nm  Nd:  YAG  lidar  with  an  optics  system  that  scans  from  horizon  to  horizon.  It  was  operated  on  1 

December,  1989  for  two  hours  when  transmissive  cirrus  were  present.  The  VIL  performed  a  pair 

of  scans  in  2  minute  cycles.  The  cycle  consisted  of  a  cross  wind  scan  (north  to  south)  followed  by 

an  along  the  wind  scan  (east  to  west).  Examples  of  images  made  from  these  scans  are  shown  in 

Figures  1 1  and  12  The  cross  wind  scans  are  the  upper  panel  and  the  along  wind  scans  are  the 
lower  panel. 

The  two  hour  period  from  19:30  to  21:30  UT  (13:30  to  15:30  CST)  began  with  light  broken 
cirrus  to  the  north  and  one  line  of  precipitating  cirrus  to  the  south  (see  Fig.  11).  The  cirrus  to  the 
north  were  in  small  patches  at  a  variety  of  levels  from  6.5  to  9  km  in  altitude.  Light  virga  were 
observed.  One  large  contrail  appeared  at  7.7  to  8.5  km  altitude,  15  km  south  of  the  VIL  (center 
point).  A  dense  line  of  precipitating  citrus  was  also  found  30  to  43  km  south  of  the  VIL.  Later 
(1.25  hours)  dense  cirrus  moved  over  the  VIL  (Fig  12).  A  nearly  continuous  cloud  from  6.5  to  93 
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km  altitude  was  present  within  most  of  the  range  of  the  VIL.  Some  optically  dense  layers  were 
found  at  7  and  9  km. 

To  quantify  the  effects  of  cloud  spatial  variability  on  scanning  method  and  FOV  size,  the 

backseat, er  measured  by  the  VIL  was  vertically  integrated  from  the  series  of  slant  angle  scans  in 

each  cross  section.  Although  the  backscatter  to  extinction  relationship  may  not  be  simple,  we  can 

assume  that  the  backscatter  is  nearly  linearly  related  to  extinction  for  these  optically  thin  clouds. 

The  vertical  integral  of  backscatter  through  the  cloud  is  proportional  to  the  optical  depth  of  the 

cloud  when  the  optical  depth  is  small.  The  effect  of  this  simplication  is  that  the  optical  depth 

structural  variations  are  most  likely  larger  than  the  backscatter  structural  variations  that  are  shown 
in  Figures  13  and  14. 

The  cross  wind  scans  were  used  to  describe  the  areal  structure  of  the  cloud.  Wind 

advection  through  each  scan  implied  that  the  clouds  were  sampled  at  approximately  3  km  intervals. 

Since  the  VlL  cross  sections  have  100  m  resolution  along  the  scan,  the  cloud  volume  was  sampled 
with  a  resolution  of  0.1  x  3  km. 

To  estimate  the  effects  of  cloud  structural  variance  on  the  averaged  vertically  integrated 

backscatter,  two  spatial  averages  were  made.  One  average  was  made  from  the  cross  wind  scans 

usmg  wind  advection  to  represent  the  satellite’s  spatial  sampling.  The  other  average  was  made 

from  the  along  wind  scans  to  represent  the  sampling  of  a  vertically  pointing  lidar  that  uses  wind 

advection  for  cloud  sampling.  Each  average  was  constructed  for  simulated  satellite  FOV 

diameters  of  1,  4,  8  and  20  km.  The  along  wind  scans  have  only  the  width  of  the  lidar  beam  ( <  1  m) 

and  thus  the  average  is  a  line  sample.  The  difference  between  these  two  averages  of  vertically 

integrated  backscatter  was  used  to  approximate  the  sampling  error  between  the  satellite  and 

stationary  lidar  sensors.  This  error  is  expressed  as  a  fraction  of  the  along  wind  lidar  average  in 

Figures  13  and  14.  To  estimate  the  possible  error  from  mis-registration  of  the  satellite  image,  the 

cross  wind  sampled  segment  was  shifted  from  the  center  position  of  the  VIL,  see  Figures  13  and 
14. 
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n.e  error  between  line  and  areal  averaged  measurements  is  the  left  most  point  m  Figures 

13  and  14.  This  error  ranged  from  5  to  29%.  It  was  higher  for  the  Ugh,  and  scatter  cirrus 

simulating  a  20  km  FOV  (Figs.  1 1  and  13)  than  the  smaller  simulated  FOVs  and  the  denser  cirrus 

(Ftgs.  12  and  14).  In  the  denser  cirrus,  the  line  vs.  area  sampling  error  ranged  from  20  to  25  % 
(Fig  14). 

The  misaUgnment  test  shows  tha,  the  area-line  error  increases  with  cross  wind 

misalignment  The  Ugh,  cirrus  (Fig.  ,3)  have  vety  large  errors  because  of  the  spacurg  between 

Cirrus  features  ,2-20  Urn  and  Cose  absent  of  the  structures  with  the  direcion  of  the  wind. 

The  error  or  dtfference  from  the  cross  wind  average  often  exceeded  the  along  wind  average  (error 

>  100%).  These  errors  were  smaUer  in  the  denser  cirrus  (Fig.  !4).  Operational  satelU.es  have 

.yptcal  nnage  registration  errors  of  12  to  20  km.  A  12  km  misaUgnment  when  compared  a 

round  based  sensor  can  cause  an  error  of  50  %  or  renter.  For  larger  FOV  sizes  of  8  and  20  km 

dus  error  decreases  sUghtly.  Because  of  this  problem,  special  effort  was  made  to  aUgn  the  sateUite 
images  to  a  one  pixel  accuracy. 


5.  DISCUSSION 


In  general,  our  observations  of  cirrus  conformed  to  the  expected  relationship  of  the  visible 
opt. cal  depth  being  twice  the  infrared  optica,  depth.  The  scatter  around  the  Wrt  .  2,  optical 

depth  ratio  mostly  came  from  the  difference  between  the  Udar  and  sateUite  methods  of  probing 
•he  clouds.  The  high  spatial  variabiUty  of  the  Coudsnnd  the  difficulty  of  evaluating  smaU  optical 
depths  from  sa, elUte  imagery  were  found  to  be  the  large,  source  of  error  in  this  comparison 
thin  hue  sampling  of  the  HSRL  vs.  the  area,  scan  of  the  sateUite  sensors  could  cause 
differences  of  20-25%  i„  measurements  of  the  same  variable.  Any  mis-regis, ration  of  the  location 
o  .  e  HSRL  tn  the  sa, elUte  image  caused  potentiaUy  larger  errors  and  had  to  be  corrected  in  this 
study.  The  initial  evaluation  of  the  HSRL  and  sateUite  data  had  more  scatter  than  shown  h,  Figure 
2.  HUS  gave  the  unpression  tha,  variations  in  cloud  structure  and  microphysical , 


i  content  were 
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causing  variations  in  the  visible-infrared  optical  depth  relationship.  Re-analysis  of  the  data  with  an 
intensive  effort  to  remove  location  uncertainties  reduced  the  scatter  in  Figure  2. 

The  satellite  measurement  of  r  a  for  optically  thin  clouds  also  had  a  bias  because  of  the 
way  in  which  the  numerator  of  Equation  6  (T4clr-T4sat)  was  evaluated.  This  required  a 
measurement  of  the  difference  between  the  clear  and  cloudy  radiances  in  the  satellite  image.  For 
optically  thin  clouds,  variations  in  the  underlying  land  surface  were  of  nearly  the  same  magnitude 
as  (Tcid'Tsat)-  The  precision  and  noise  characteristics  of  the  satellite  sensors  require  that  the 
numerator  terms  in  Equation  6  be  at  least  0.8  K  different  (T^-T^  >  0.8K)  for  the  cloud  to  be 
distinguished  on  the  satellite  image.  This  forces  a  minimum  of  0.03-0.05  for  tu.  Any  error  in  the 
terms  in  Equation  6  from  spatial  variance  tend  to  increase  r  a  because  errors  that  would  reduce  tu 
make  cloud  detection  impossible.  This  is  possibly  why  7  of  the  9  points  in  Figure  2  with  <  1 

appear  above  the  2.0  line.  It  also  indicates  that  the  minimum  optical  depth  (infrared)  detectable  in 

y 

,,  a  satellite  image  is  about  0.05. 

Optically  thin  cirrus  (t™  <  1)  best  conformed  to  the  assumption  made  by  most  satellite 
cirrus  cloud  detection  algorithms  that  cirrus  clouds  can  be  described  as  a  single  thin  layer  even 
though  they  span  depths  of  1-3  km.  This  assumption  is  possible  when  the  clouds  have  little  vertical 
variance  in  radiative  properties. 

For  dense  cirrus  with  optical  depths  >  1  (t„  >  0.5),  the  cirrus  were  often  over  3  km 
thick  and  in  two  or  more  layers.  Vertical  variations  in  optical  density  were  common  and  liquid 
water  layers  were  often  present  inside  the  ice  clouds.  The  retrieval  of  a  single  cloud  layer  from  the 
H3RS  data  to  describe  the  infrared  characteristics  of  these  clouds  appears  to  work  as  a  first 
approximation  of  infrared  divergence.  More  detail  in  the  vertical  would  be  helpful  for  the 
development  of  radiative  transfer  models.  But  the  fact  that  the  infrared  optical  depth  calculated 
from  the  single  cloud  level  solution  followed  the  theoretical  relationship  of  being  nearly  one  half  of 
the  visible  optical  depth  is  very  useful  for  inter-relating  cloud  climatologies.  Development  of 
radiative  transfer  algorithms  for  global  numerical  models  requires  global  data  sets  on  clouds  and 
their  radiative  properties.  Satellites  are  currently  the  best  source  of  this  data. 
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We  found  super  cooled  liquid  water  layers  in  32%  of  the  48  cirrus  cases  measured  by  the 
HSRL  which  was  more  frequent  than  we  expected.  These  case  include  an  extra  25  days  which 
were  not  analyzed  with  satellite  data.  This  data  set  may  not  be  large  enough  to  state  the  frequency 
of  liquid  water  in  all  cirrus  clouds.  The  HSRL  was  operated  mostly  during  NOAA  12  overpasses  at 
dusk  and  after  dusk.  Cirrus  data  could  only  be  obtained  when  obstructions  from  lower  water 
clouds  were  not  present  or  minimal  enough  to  permit  the  view  of  the  cirrus.  Cirrus  (ice  clouds) 
obviously  exist  above  low  altitude  liquid  water  clouds  but  we  have  little  data  from  these  cirrus 
because  of  low  cloud  obscurations. 

6.  CONCLUSIONS 

The  assumption  of  a  single  cloud  level  used  by  the  satellite  cloud  climatologies  is  applicable 
to  a  largejortion  of  the  satellite  data.  The  manner  in  which  optical  depths  are  calculated  from  the 
10:8  Mm  window  channel  measurements  appears  to  approximate  the  vertically  averaged  optical 
depths  for  a  variety  of  cloud  structures  which  are  far  more  complicated  than  one  layer. 

The  cirrus  climatologies  derived  from  different  satellite  data,  visible  and  infrared  10.8  nm 
window,  can  be  related  using  the  2.0  optical  depth  ratio  to  form  longer  records  of  cirrus  cloud 
cover  and  describe  inter-annual  changes.  This  work  is  needed  to  asses  the  relationship  of  cirrus 
clouds  to  climatic  changes  and  their  role  in  the  earth’s  heat  budget.  This  study  used  clouds  of 
visible  optical  depths  ranging  from  0.25  to  3.0.  According  to  Wylie  et  al.  (1994),  these  clouds  cover 
32%  of  the  earth  which  is  a  significant  part  of  the  earth’s  cloud  cover. 

Super  cooled  liquid  water  layers  are  common  inside  ice  clouds  and  can  easily  change  the 
radiative  characteristics  of  these  clouds.  They  were  found  in  nearly  one  third  of  our  cirrus  data. 
Their  impact  on  radiative  processes  needs  to  be  evaluated  and  models  of  radiative  processes  in 
clouds  need  to  include  them. 

The  problem  of  matching  point  or  line  sensors  (lidars)  to  areal  scanners  (satellites)  is  an 
obstacle  to  field  programs  and  the  use  of  their  data  with  global  cloud  climatologies.  Miss-matches 
in  sampling  and  scanning  can  lead  to  large  errors  when  data  sets  are  compared  and  these  errors 
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can  be  mistaken  for  variances  in  radiative  properties  of  the  clouds.  Scanning  sensors  and  a  large 
number  of  cases  are  needed  from  these  field  programs  to  define  cloud  radiative  properties. 
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TABLE  1.  The  dates  and  times  of  the  satellite-lidar  comparisons. 


Month/Date/Year 

Time 

Satellite 

08/24/93 

01:30  UTC 

NOAA  12 

09/02/93 

01:37 

NOAA  12 

09/08/93 

01:08 

NOAA  12 

09/09/93 

00:46 

NOAA  12 

10/02/93 

00:51 

NOAA  12 

10/12/93 

00:36 

NOAA  12 

10/25/93 

00:57 

NOAA  12 

10/26/93 

00:35 

NOAA  12 

10/28/93 

01:32 

NOAA  12 

11/10/93 

00:12 

NOAA  12 

11/11/93 

01:31 

NOAA  12 

02/04/94 

00:06,  00:36, 

01:06  and  01:36 

GOES 

02/15/94 

01:06 

GOES 

02/17/94 

01:06 

GOES 

03/04/94 

01:06 

GOES 

03/08/94 

01:09 

NOAA  12 

03/17/94 

01:09 

GOES 

03/26/94 

00:06,  00:36, 

GOES 

04/07/94  > 

01:06  and  01:36 

00:22 

NOAA  12 

,04/11/94 

00:36 

NOAA  12 

04/15/94 

00:49 

NOAA  12 

Multiple  layers 
Thin  and  variable 
Multiple  layers 
Thick  cirrus 
Thick  cirrus 
Thin 

Thin  and  variable 
Uniform  layer 
One  thick  layer 
Two  layers 
Two  layers 
One  thick  layer 

Thin  layer 
Thin  layer 
Thin  layer 
Two  layers 
Thin  layer 
Thick  layer 

Thin  layer 
Thick  layer 
Two  layers 
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Figure  1:  The  physical  thickness  vs.  the  visible  optical  depth  of  semi-transparent  clouds  observed 
by  the  HSRL  in  1993  and  1994. 


Figure  2:  Scatter  plot  of  the  visible  and  infrared  optical  depths  measured  by  the  combination  of  the 
HSRL  lidar  (visible  optical  depth  at  0.53  /xm)  and  NOAA  and  GOES  satellite  data  (infrared  10.8 
pm  window  channel).  The  diagonal  line  is  a  ratio  of  2.0  (visible/IR).  Error  bars  are  shown  on 
each  data  point. 

Figure  3:  Time-height  image  of  the  backscatter  cross  section  pa(x)  ^(seeEq.2) 

from  the  HSRL  (upper  panel).  The  lower  panel  is  the  visible  optical  depth  integrated  from  8.5  to 
12  km  altitude  on  26  October  1993. 


Figure  4:  NOAA  12  1  km  HRPT  image  from  AVHRR  Channel  4  (10.8  pm)  from  26  October  1993. 
This  image  was  remapped  into  an  equal  distant  latitude-longitude  projection.  The  location  of  the 
HSRL  is  indicated  by  the  black  dot. 


Figure  5:  Time-height  image  of  the  backscatter  cross  section  pa(r )  (see  Eq.  2) 

from  the  HSRL  (upper  panel).  The  lower  panel  is  the  visible  optical  depth  integrated  from  4  to  10 
-  km  altitude  on  11  November,  1993.  The  extremely  white  layers  between  5  and  5.5  km  were  liquid 
water  as  determined  from  depolarization  measurements  (not  shown). 


Figure  6:  NOAA  12  1  km  HRPT  image  from  AVHRR  Channel  4  (10.8  jum)  from  11  November 
1993.  This  image  was  remapped  into  an  equal  distant  latitude-longitude  projection.  The  location 
of  the  HSRL  is  indicated  by  the  black  dot. 


Figure  7:  Time-height  image  of  the  backscatter  cross  section  pa{ r)  (see  Eq.  2) 

from  the  HSRL  (upper  panel).  The  lower  panel  is  the  visible  optical  depth  integrated  from  7  to  11 
km  altitude  on  2  September,  1993. 


Figure  8:  NOAA  12  AVHRR  Channel  4  infrared  image  from  2  September,  1993. 

Figure  9:  Time-height  image  of  the  backscatter  cross  section  y9fl(r)  (see  Eq.  2) 

from  the  HSRL  (upper  panel).  The  lower  panel  is  the  visible  optical  depth  integrated  from  3.5  to 
11  km  altitude  on  15  April,  1994. 


Figure  10:  NOAA  12  AVHRR  Channel  4  infrared  image  from  15  April,  1994. 

Figure  11:  The  cross  wind  (top  panel)  and  along  the  wind  (bottom  panel)  images  of  range 
corrected  backscatter  from  the  Volume  Imaging  Lidar  for  19:50  UT  (13:50  CST)  on  1  December 
1989  at  Madison,  WI. 
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Figure  12:  Same  as  Fig.  11  for  20:07  UT  (15:07  CST)  o„  1  December  1989. 

volumes  for  the  period  of  19:30  to  20:30  UT,  1  December,  1989. 

Figure  14:  Same  as  Fig.  12  for  the  period,  20:30  to  21:30  UT,  1  December,  1989. 
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Figure  1:  The  physical  thickness  vs.  the  visible  optical  depth  of  semi-transparent  clouds  observed 
by  the  HSRL  in  1993  and  1994. 
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Figure  3:  Time-height  cross  section  of  backscatter  (upper  panel)  and  visible  optical  depths  between 
8.5  and  12  km  (lower  panel)  from  the  HSRL  on  26  October  1993. 


Figure  4:  NOAA  12  1  km  HRPT  image  from  AVHRR  Channel  4  (10.8  /um)  from  26  October  1993. 
This  image  was  remapped  into  an  equal  distant  latitude-longitude  projection.  The  location  of  the 
HSRL  is  indicated  by  the  black  dot. 
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Figure  6:  NOAA 12  1  km  HRPT  image  from  AVHRR  Channel  4  (10.8  /un )  from  11  November 
1993.  This  image  was  remapped  into  an  equal  distant  latitude-longitude  projection.  The  location 
of  the  HSRL  is  indicated  by  the  black  dot. 


Figure  6:  NOAA  12  1  km  HRPT  image  from  AVHRR  Channel  4  (10.8  /im)  from  11  November 
1993.  This  image  was  remapped  into  an  equal  distant  latitude-longitude  projection.  The  location 
of  the  HSRL  is  indicated  by  the  black  dot. 
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Figure  7:  Time-height  cross  section  of  backscatter  (upper  panel)  and  visible  optical  depth  between 
7  and  11  km  (lower  panel)  from  2  September,  1993,  measured  by  the  HSRL. 


Figure  8:  NOAA  12  AVHRR  Channel  4  infrared  image  from  2  September,  1993. 
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^9V1T^height  crosfNs5ction  of  backscatter  (upper  panel)  and  visible  optical  depth  between 
3.5  and  1 1  km  (lower  panel)  from  15  April,  1994,  measured  by  the  HSRL. 


Figure  10.  NOAA  12  AVHRR  Channel  4  infrared  image  from  15  April,  1994. 
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Figure  14:  Same  as  Fig.  12  for  the  period,  20:30  ,o  21:30  UT,  1  December,  1989. 


